AD-A1 62  119  OPTIMIZATION  OF  TIP  5T0RE  MODEL ING<U>  OKLAHOMA  UNIV 
NORMAN  SCHOOL  OF  AEROSPACE  MECHANICAL  AND  NUCLEAR 
ENGINEERING  A  G  STRIZ  ET  AL  MAR  85  AFOSR-TR-85-1889 
UNCLASSIFIED  AFOSR-83-8184  F/G  28/4 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


AFOSR-TR- 


p 


0  39 


SECURITY  CLASSIFICATION  OF  THIS  PAGt 


la  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


3a.  SECUAITV  CLASSIFICATION  AUTHORITY 


This  PAGE  l 


-  J  i  4  :• 

.  I  *>  ^4,-*  '  * 


:  c. 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3.  OlSTRIBUTION/AVAILABH 


3b  OECLASSIFICATION/OOINNGRADING  SCHEDULE 


Approved  fer  public  release^ 
dlstributica  ualimitBd 


«.  PERFORMING  ORGANIZATION  REPORT  NUMBERISI 

5.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 

AFOSRTR- 

6a  NAME  OF  PERFORMING  ORGANIZATION 

UNIVERSITY  OF  OKLAHOMA 

Sb.  OFFICE  SYMBOL 
< If  applicable  I 

7«.  NAME  OF  MONITORING  ORGANIZATION 

AFOSR/NA 

6c.  ADDRESS  (City.  State  and  ZIP  Coda 

SCHOOL  OF  AEROSPACE,  MECHANICAL 
NORMAN,  OKLAHOMA  73019 

&  NUCLEAR  ENGR 

?b.  AOORESS  (City.  State  and  ZIP  Code’ 

BOLLING  AFB  DC  20332-6448 

Ba  name  of  funding/sponsoring 
organizationa1r  FORCE 

OFFICE  OF  SCIENTIFIC  RESEARCH 

8b.  OFFICE  SYMBOL 
(If  applicable) 

AFOSR/NA 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

AFOSR-83-0184 

Be.  AOORESS  (City.  State  mid  ZIP  Code) 

BOLLING  AFB,  DC  20332-6448 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM  PROJECT 

ELEMENT  NO.  NO. 


11.  TITLE  (Include  Security  Claeti  fleet  ion) 

OPTIMIZATION  OF  TIP  STORE  MODELING  (UNCLASSIFIED) 


12.  personal  authorisi 

ALFRED  G  STRIZ.  SUNG-KUK  JANG 


ISA  type  of  REPORT 


61102F 


WORK  UNIT 
NO 


1 13b.  TIME  COVERED 

1*.  DATE  OF  REPORT  fYr..  Mo..  Day) 

1 

1985  March 

COSAT  I  CODES 


GROUP  I  SUB.  GR. 


18.  SUBJECT  TERMS  fConWnue  on  reverie  if  neeeuary  and  identify  by  block  number i 

TIP  STORE  .WING,  MODELING,  UNSTEADY  AERODYNAMICS 


ABSTRACT  /Continue  on  reverie  if  neceuar y  and  identify  by  block  number i 

~s  Due  to  the  high  cost  of  performing  detailed  three-dimensional  flutter  analyses  for  air¬ 
craft  which  are  carrying  large  numbers  and  types  of  external  stores,  it  is  often  not  econom¬ 
ically  feasible  to  include  sophisticated  store  aerodynamics  especially  when  they  will  not 
cause  significant  changes  to  the  flutter  results.  Thus,  the  stores  are  presently  often 
modeled  as  simple  flat  plates.  This  tends  to  decrease  the  accuracy  of  the  computed  pressure 
distributions  as  the  spanwise  flow  field  and  the  wing-body  interaction  are  not  represented 
correctly. 

In  this  report,  therefore,  two  methods,  a  doublet  lattice  method  for  subsonic  flow^ 
(computer  code  H7WC  by  Giesing,  Kalman,  and  RoddenV^and  a  kernel  function  method  for  subsonic 
and  supersonic  flow  ^(computer  code  ANKF  by  Cunningham)^are  used  for  the  computation  of  un¬ 
steady  pressure  distributions  and  forces  on  an  F-5  wing  with  a  tip  mounted  launcher/store 
combination  for  various  store  models  consisting  of  simple  flat  plates,  endplated,  and 
axisymmetric  bodies.  The  resulting  aerodynamic  data  are  compared  to  existing  experimental,,  . 


30.  OISTRI BUTION/AVAI LABILITY  OF  ABSTRACT 
UNCLASSIFICO/UNLIMITEO  0  SAME  AS  RPT.  □  OTIC  USERS  □ 


33a  NAME  OP  RESPONSIBLE  INC  VIOUAL 

S  D  WILSON 


31.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


32b  TELEPHONE  NUMBER  22c.  OFFICE  SYMBOL 

I Include  Area  Code) 

767-4935  AFOSR/NA 


jjCUWlTy  classification  op  THIS  page 


->.nd  numerical  results  for  the  sate  «ing/store  combination. 

,  dyne„I2C£”«:  fmlicit,  end  accuracy  with  respect  to  the  eero- 

increase  in  complexity  i  e  nanol  n  k  °  y  n  ®ract*on  was  improved  without  a  considerable 
n  complexity,  i.e.  panel  number,  and  thus  computational  effort.  fS 

.  /  .  >  *  '  1  *  '  a  -  -d-V-.d  . 


o  d 


/  f 


C;/> 


OPTIMIZATION  OF  TIP  STORE  MODELING 


Alfred  G.  Striz  and  Sung-Kuk  Jang 

University  of  Oklahoma 
Norman,  OK  73019 


March  1985 


Technical  Report 

Final  Report  for  Period  May  1983  -  September  1984 


Approved  for  public  release;  distribution  unlimited. 


Air  Force  Office  of  Scientific  Research 

Air  Force  Systems  Command 

Bolling  Air  Force  Base,  D.C.  20332 


FOBEWOBD 


This  report  was  prepared  by  Dr.  Alfred  G.  Striz  of  the 
School  of  Aerospace,  Mechanical  and  Nuclear  Engineering  at  the 
Oniversity  of  Oklahoma  in  Norman,  OK,  under  AFOSE  Grant  83-0184, 
"Optimization  of  Tip  Store  Modeling".  The  research  was  adminis¬ 
tered  by  Captain  Dwight  R.  McGhee  and  the  program  manager  was 
Dr.  James  D.  Wilson,  both  of  the  Air  Force  Office  of  Scientific 
Research  at  Bolling  Air  Force  Base,  D.C. 


The  report  covers  work  conducted  during  the  grant  period 
from  May  1983  to  September  1984.  Dr.  Alfred  G.  Striz  was  the 
principal  investigator,  Sung-Kuk  Jang  was  the  graduate  research 
assistant.  Bob  Bunton  and  Dave  Eelk  of  the  Air  Force  Armament 
Laboratory  at  Eglin  Air  Force  Base,  FL,  and  Herbert  Cunningham 
of  the  NASA-Langley  Research  Center  in  Hampton,  VA,  were  helpful 
in  providing  needed  information. 

Sung-Kuk  Jang  will  receive  his  M.S.  degree  for  this  re¬ 
search  in  March  of  198b  with  a  thesis  entitled  "Optimization  of  j 

Wing  Tip  Store  Modeling  for  Aircraft  in  Subsonic  and  Supersonic 

Flow".  He  presented  partial  results  of  this  work  previously  at 

the  1984  meeting  of  the  Oklahoma  Academy  of  Sciences  in  Ada,  OK, 

-  iii  - 


'■  c- 

■  -O  1.1 


‘  .'-Vi  . 


and  at  the  1985  AI AA/ ASME  Symposium  VII  in  Norman,  OK. 

A  paper  with  the  title  "Optimization  of  Tip  Store  Mode¬ 
ling"  by  Dr.  Alfred  G.  Striz  and  Sung-Kuk  Jang  will  be  submitted 
to  AIAA  for  possible  publication  in  the  Journal  of  Aircraft. 


TABLE  OF  CONTENTS 


FOREWORD  .  .  -  . 
TABLE  OF  CONTENTS 
LIST  OF  SYMBOLS  . 
LIST  OF  FIGURES  . 
LIST  OF  TABLES  . 


Chapter 

I.  INTRODUCTION  . 

Background  . 

Scope  of  Present  Research 

II-  DEVELOPMENT  OF  THEORY  .  .  .  . 

Doublet  Lattice  Method  .  . 
Introduction  .  ...  . 
Theoretical  Development 
Kernel  Function  Method  .  . 
Introduction  .  .  .  .  . 
Theoretical  Development 

III-  NUMERICAL  MODELING  . 

Doublet  Lattice  Method  .  . 
Kernel  Function  Method  .  . 

IY.  RESULTS  . 

Clean  Wing  Case  . 

Doublet  Lattice  Method  .  . 
Kernel  Function  Method  .  . 
Comparison  ........ 

Y.  CONCLUSION  AND  FUTURE  WORK  . 

Conclusion  ..  . 

Future  Work  . . 


REFERENCES 


LIST  OF  SYMBOLS 


Speed  of  sound 

Coefficient  in  the  approximation  of  the 
vibration  inodes 

Wing  semi-chord  nondiaensionalized  by  C 
Pressure  coefficient 
Boot  chord 

Normal  force  coefficient 
Kernel 

ti>C  n 

Reduced  frequency  (*  —  ) 

2Um 

Coordinates  of  control  points 
Mach  number 
Lifting  pressure 
weighting  functions 
Dynamic  pressure  (*  ^pU^2  ) 

Wing  semi-span 

Free-stream  velocity 

Flow  velocity  vector  (  =V(x,y,z,t)  ) 

Onnormalized  normalvash 

W/(J 

oo 

Modal  displacement 

Coordinates  of  receiving  points 

Coordinates  of  load  (integration)  point; 


Coordinates  of  the  load  points  in  the 


transformed  plane 
Lifting  pressure  coefficient 
Coordinates  of  sending  points 

Lateral  coordinate  in  the  plane  of  the  surface 
Density  of  air  at  sea  level 
Velocity  potential  [  =Q  (x,y,z,t)  ) 


Frequency  of  oscillation 
Sweeo  angle 


LIST  OP  FIG0H2S 


Piqure  page 

1.  Configuration  of  F-5  Ming  .............10 

2.  Configuration  of  P-5  Ming  with  Stores  .......  11 

3.  Ring  and  Element  Coordinates . 13 

4.  Idealization  of  Rihg  Panel  into  Boxes  .......15 

5.  Coordinates  of  Transformed  Surface  .........20 

6.  Surface  and  Loading  Types . 23 

7.  Various  Tip  Store  Cross  Sections  ....  25 

8.  Mode  Shape  of  the  Clean  Ring  (3=0. y5,  F=20Hz,  a  =0)  .  27 

9.  Geometrical  Model  of  Ring-Launcher  (DIM) . 31 

10.  Geometrical  Model  of  Missile  (DLM)  . . 32 

11.  Geometrical  Model  of  Ring- Launcher  (KFH)  ......  33 

12.  Geometrical  Model  of  Missile  (KFM)  ..........34 

13.  Unsteady  Lift  Coefficient  Distributions  of  Clean 

Ring  . . 38 

14.  Unsteady  Lift  Coefficient  Distributions  (3=0.6, 

k=0.200,  DLM ,  22  Models) . 40 

15.  Unsteady  Lift  Coefficient  Distributions  (3=0.7, 

k=0. 173,  DLM,  22  Models) . 41 

16.  Unsteady  Lift  Coefficient  Distributions  (M=0.8, 

k=Q .  1 5 3 ,  DLM,  22  Models) . 42 

17.  Unsteady  Lift  Coefficient  Distributions  (M=0.9, 

)c=0. 138,  DLM,  22  Models) . 43 

18.  Unsteady  Lift  Coefficient  Distributions  (3=0.6, 

k=0.200,  DLM,  8  Models) . 4b 


v  i  f  i 


Unsteady  Lift  Coefficient  Distributions  (M=0.7, 
k=0.173,  DLM,  8  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (M=0.8, 
*=0.153,  DLM,  8  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (M=0.9, 
*=0.138,  DLM,  3  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (M=0.6, 
*=0.200,  KFM,  6  Models) . 

Unsteady  Lift  Coefficient  Distributions  (M=0.7, 
*=0.173,  KFM  ,  b  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (M=0.8, 
*=0.153,  KFM,  6  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (8=0.9, 
*=0.138,  KFM,  6  Models)  . 

Unsteady  Lift  Coefficient  Distributions  (3=1.1, 
*=0.116,  KFM,  6  Models)  .  -  . 

Unsteady  Lift  Coefficient  Distributions  (M=i.3b, 
*=0.101,  KFM,  6  Models)  . 

Pressure  Distributions  at  Section  2  (Beal,  M=0.8 
*=0.153)  for  Various  Models  ..  . 

Pressure  Distributions  at  Section  2  (Iaag,  M=0.8 
*=0.153)  for  Various  Models  .  . 

Pressure  Distributions  at  Section  5  (Beal,  M=0.8 
*=0.153)  for  Various  Models  ......... 

Pressure  Distributions  at  Section  5  (Iaag,  H=0.8 
*=0.153)  for  Various  Models  ......... 

Pressure  Distributions  at  Section  8  (Beal,  H=0.8 
*=0.153)  for  Various  Models  .  .... 

Pressure  Distributions  at  Section  8  (Iaag,  N=0.8 
*=0.153)  for  Various  Models  . 


LIST  OP  TABLES 


Table  page 


1.  Coefficients  "amn"  for  Par ions  Vibration  Modes  ...  37 

2.  Onsteady  Lift  Coefficient  Distributions  of  Clean 

Sing  ......................37 

3.  Onsteady  Lift  Coefficient  Distributions  (Heal, 

3=0.6,  *=0 . 200 ,  DL3 ,  8  Models) . 45 

4.  Onsteady  Lift  Coefficient  Distributions  (Imag, 

3=0.6,  X=0 . 200 ,  DL3,  8  Models)  .  45 

5.  Onsteady  Lift  Coefficient  Distributions  (Real, 

3=0.7,  X=0. 173,  DL3,  8  Models)  .  48 

6.  Onsteady  Lift  Coefficient  Distributions  (Imag, 

3=0.7,  X=0 . 1 73,  DLS ,  8  Models) . 48 

7.  Onsteady  Lift  Coefficient  Distributions  (Real, 

3=0.8,  *=0.153,  DL3 ,  8  Models)  .  50 

8.  Onsteady  Lift  Coefficient  Distributions  (Imag, 

3=0.8,  X=0. 153,  DLfl,  8  Models)  .  50 

9.  Onsteady  Lift  Coefficient  Distributions  (Heal, 

3=0.9,  k=0 . 138 ,  DL3,  8  Models) . 52 

10.  Onsteady  Lift  Coefficient  Distributions  (Iaag, 

3=0.9,  *=0.138,  DL3,  8  Models)  .  52 


11.  Comparison  of  Effort  (3=0.6,  *=0.200,  DL3 ,  8  Models)  55 

12.  Comparison  of  Effort  (M=0.7,  *=0.173,  DLM ,  8  Models)  55 

13.  Comparison  of  Effort  (3=0.8,  *=0.153,  DLM,  8  Models)  5b 

14.  Comparison  of  Effort  (3=0.9,  X=0.138,  DLM,  8  Models)  5b 


15.  Onsteady  Lift  Coefficient  Distributions  (Real, 

3=0.6,  *=0.200,  K?3  ,  6  Models) . 58 

16.  Onsteady  Lift  Coefficient  Distributions  (Imag, 

3=0.6,  *=0.200,  KFH ,  6  Models) . 58 


17.  Unsteady  Lift  Coefficient  Distributions  (Heal, 

8=0.7 ,  *=0 .173,  KFM,  6  Models) . 60 

18.  Unsteady  Lift  Coefficient  Distributions  (Imag, 

8=0.7,  *=0. 173,  KFM,  6  Models)  .  60 

19.  Unsteady  Lift  Coefficient  Distributions  (Real, 

M=0.8,  *=0.153,  KFM,  6  Models)  .  62 

20.  Unsteady  Lift  Coefficient  Distributions  (Iaag, 

8=0.8,  *=0.153,  KFM ,  6  Models) . 62 

21.  Unsteady  Lift  Coefficient  Distributions  (Heal, 

M=0.9,  *=0.138,  KFM,  6  Models)  .  64 

22.  Unsteady  Lift  Coefficient  Distributions  (Imag, 

8=0.9,  *=0.138,  KFM,  6  Models)  .  64 

23.  Unsteady  Lift  Coefficient  Distributions  (Beal, 

M=1.1,  *=0.116,  KFM,  6  Models)  .  . . 67 

24.  Unsteady  Lift  Coefficient  Distributions  (Imag, 

.1  =  1.1,  *=0.116,  KFM,  6  Models) . 67 

25.  Unsteady  Lift  Coefficient  Distributions  (Heal, 

1=1 .35 ,  *=0.101,  KFM,  6  Models) . 69 

26.  Unsteady  Lift  Coefficient  Distributions  (Imag, 

1=1 .35 ,  *=0.101,  KFM,  6  Models) . 69 


27.  Comparison  of  Effort  (3=0.6,  *=0.200,  KFM,  6  Models)  72 

28.  Comparison  of  Effort  (H=0.7,  k=0.173,  KFM,  6  Models)  72 

29.  Comparison  of  Effort  (M=0. 8,  *=0.153,  KFM,  6  Models)  73 

30.  Comparison  of  Effort  (M=0.';,  *=0.138,  KFM,  6  Models)  73 


31.  Comparison  of  Effort  (8=1.10,  *=0.116,  KFM,  6 

Models)  . . 74 

32.  Comparison  of  Effort  (M=1.35,  *=0.101,  KFM,  6 

Models)  ..................  ...74 

33.  Pressure  Distributions  at  Section  2  (Real,  M=0.8, 

*=0.153)  for  Various  Models  . . 75 

34.  Pressure  Distributions  at  Section  2  (Imag,  3=0.8, 

k=0.l53)  for  Various  Models  .  75 

35.  Pressure  Distributions  at  Section  5  (Real,  8=0.8, 

*=0.153)  for  Various  Models . 78 


xt  - 


cl 

y 

,v 

■ 


■ 


;v 

.V 


Chapter  I 
IHTBODUCTIO N 

When  large  external  bodies  or  stores,  such  as  missives, 
bombs,  or  tael  tanks,  are  added  to  the  wings  of  an  airplane, 
the  dynamic  characteristics  of  the  aircraft  will  be  changed. 
In  particular,  the  flutter  speed  of  the  aircraft  may  be 
significantly  affected  because  of  the  inertial,  elastic,  and 
aerodynamic  interaction  between  the  wing  and  its  mounted 
stores.  Thus,  theoretical  aeroelasticity  is  essential  for 
the  flutter  evaluation  of  this  type  of  airplanes.  ' 

However,  until  recently,  the  computation  of  aerodynamic 
loadings  over  complex  geometric  configurations  was  not 
possible.  Fortunately,  the  development  of  three-dimensional 
aerodynamic  codes  capable  of  calculating  wing-body 
interaction  has  now  advanced  to  the  point  where  the  flow 
about  these  complex  structures  and  thus  the  aerodynamic 
forces  and  pressures  can  be  predicted  quite  accurately.  In 
order  to  achieve  this  accuracy,  rather  detailed  models  nave 
to  be  established. 
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Several  recent  studies  and  tests  have  considered  the 
effects  of  store  aerodynamics  on  unsteady  wing  airloads  with 
emphasis  placed  on  wings  with  tip  missiles. 

Chadwick  (1)  presented  a  nonplanar  lifting  surface 
theory  to  predict  sideioads  on  underwing  pylon  store  models. 
The  wing  pylon  model  was  limited  to  incompressible  flow  and 
slender  body  theory  was  used  to  model  the  store.  Kraus  (2) 
used  panel  methods  to  compute  forces  and  moments  on  external 
stores  mounted  underneath  an  F-4  wing  in  subsonic  flow. 
Plartin,  Saunders,  and  Smith  (3)  reported  an  analysis  of  the 
aerodynamic  interference  problem  for  external  aircraft 
stores  by  use  of  the  image  system  technique.  The  stores 
were  assumed  to  be  slender,  axisymmetric  bodies  and  the 
interference  was  analyzed  by  first  assuming  a  cross-flow 
solution.  In  the  Netherlands,  Eoos,  BenneJcers,  and  Zwaan 
(4)  of  NLB  described  a  panel  method  for  the  calculation  of 
aerodynamic  loadings  on  harmonically  oscillating  wing/oody 
configurations  xn  subsonic  flow.  The  loads  on  the  body  were 
represented  by  an  unsteady  source  panel  distribution,  while 
the  unsteady  aerodynamic  forces  on  the  lifting  surfaces  were 
computed  by  the  doublet  lattice  method.  A  wing/tip  tank  and 
a  pylon/store  configuration  were  evaluated.  Parker  (5)  used 
the  doublet  lattice  method  together  with  the  method  of 
images  to  predict  unsteady  aerodynamic  coefficients  and 
pressures  for  an  F-b  wing  with  various  store  configurations 
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oscillating  in  subsonic  flow.  A  variety  of  finite  element 
lattice  configurations  was  evaluated  tor  numerical 
convergence.  CenJco,  Tinoco,  Dyer,  and  DeJongh  (6)  analysed 
the  complex  flow  associated  with  aircratt/weapons  carriage 
and  mutual  interference  during  store  separation  by  means  of 
the  PAN  AIR  code.  This  program  is  based  on  a  higher  order 
panel  method.  Dusto  (7)  also  presented  a  higher  order 
paneling  method  to  compute  steady  and  unsteady  flows  aaout 
an  P-5  wing  with  and  without  underwing  and  wingtip  store. 
These  results  and  similar  ones  from  another  higher  order 
panel  method  and  a  doublet  lattice  method  were  compared  by 
Sotomayer,  Dusto,  Epton,  and  Johnson  (8)  to  experimental 
results  for  the  same  configurations  from  NLR  (9) .  In  a 
wing/store  flutter  suppression  investigation.  Noil, 
Huttsell,  and  Cooley  (10)  used  the  subsonic  doublet  lattice 
option  in  the  automated  flutter  analysis  module  of  tne 
PASTOP  program  (11)  to  perform  unaugmented  flutter 
calculations,  and  to  obtain  unsteady  aerodynamic  forces  for 
a  YF-17  wing  with  an  AIM-7S  missile  mounted  on  an  outboard 
wing  pylon  and  with  an  empty  launcher  rail.  Pollock., 
Sotomayer,  Huttsell,  and  Cooley  (12)  evaluated  various 
methods  for  the  predictions  and  prevention  of  wing/store 
flutter.  The  measured  aerodynamic  data  of  the  NLR 
experiments  were  compared  to  unsteady  aerodynamic  data  from 
a  paneling  method  by  Woodward  (13).  Flutter  analyses  were 


performed  using  the  doublet  lattice  module  in  the  flutter 
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code  FACES  (14)  for  the  same  configurations.  Triplett  (15) 

|  used  doublet  lattice  theory  to  evaluate  the  effect  of 

detailed  aerodynamic  modeling  of  the  F/A-18  wing  with 
underwing  stores  and  tip  missile.  He  presented  results 
j  showing  the  effect  of  individual  system  components  on 

flutter.  Finally,  Turner  (16)  presented  an  analytical  study 
of  the  effect  of  store  aerodynamics  on  wing/store  flutter  by 
j  using  multivariate  analysis  techniques. 

In  all  this  above  mentioned  research,  the  pylon/store 
or  launcher/store  configurations  were  modeled  in 
I  three-dimensional  detail.  For  an  applied  flutter  analysis 

in  the  flutter  certification  of  aircraft  with  stores,  this 
becomes  very  costly,  as  pointed  out  by  Turner. 

! 

1.2  SCOPE  OF  PBESEBT  BESEABCH 

The  computations  of  the  unsteady  aerodynamic 
coefficients  for  an  aircraft  are  the  single  most  costly  item 
in  any  flutter  analysis.  Also,  owing  to  the  variety  of 
stores  that  can  be  carried  on  typical  modern  tactical 
military  aircraft,  the  total  number  of  possible 
aircraft/store  configurations  is  rather  large.  Thus,  with 
thousands  of  flutter  computations  to  be  performed,  the  cost 
becomes  prohibitive. 

In  many  cases,  this  increased  cost  is  not  justified 
because  store  aerodynamics  tend  to  have  a  small  effect  on 
flutter  speed  as  shown  by  Turner  (lb).  However,  there  are 
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important  instances  whore  neglecting  store  aerodynamics  wiaa 
lead  to  a  dangerous  over-estimation  of  flutter  speed. 
Determination  when  store  aerodynamics  snould  bo  included  in 
the  flutter  analysis  is  a  major  problem. 

As  a  compromise  it  has  become  common  practice  to  model 
the  stores  solely  as  flat  plates  when  computing  the 
aerodynamic  loads  on  oscillating  wing-body  combinations. 
This,  however,  tends  to  reduce  the  accuracy  of  the  results. 

According  to  Chapman  M7),  for  panel  methods'  the 
computation  time  varies  with  the  number  of  panels  n  as 
somewhere  between  n 2  and  a3.  Thus,  a  compromise  has  to  be 
sought  between  the  complex  geometric  modeling  in  a  one-time 
evaluation  of  a  computer  code  and  the  requirement  for 
simplicity  of  the  model  to  save  cost  ^n  repeated  flutter 
analyses  of  multiple  aircraft/store  configurations,  while 
also  preserving  sufficient  accuracy. 

In  this  investigation,  studies  are  conducted  to 
evaluate  various  geometric  models  for  pylons  and  stores  such 
as  endplates,  endplate  -  flat  plate  combinations,  and  these 
models  at  various  dihedral  angles,  in  comparison  to  the 
coaputationai  flat  plate  models  now  in  use.  This  is  done 
with  respect  to  complexity  and  accuracy  using  a  kernel 
function  code,  ANKF,  by  Cunningham  ( i 8)  and  a  doublet 
lattice  code,  H7WC,  by  Giesing,  Kalman,  and  Rodden  M 9)  to 
investigate  whether  the  accuracy  of  the  aerodynamic  resuits 
can  be  increased  by  using  these  somewhat  more  complex  modeis 


Chapter  II 


DEVELOPMENT  OF  THEOBI 


Steady  and  unsteady  aerodynamic  loads  on  airplanes  are 
used,  in  general,  to  analyze  flutter,  gust,  and  freguancy 
response  as  well  as  for  the  computation  of  static  and 
dynamic  stability  derivatives.  Thus,  the  accurate  and 
efficient  prediction  of  these  aerodynamic  loads  is  an 
essential  requirement  for  any  such  analysis. 

The  use  of  lifting  surface  theory  has  fcaen  the  classic 
way  to  solve  the  unsteady  compressible  flow  problem. 
Starting  from  conservation  of  mass,  momentum,  and  energy  as 
well  as  from  Gibb's  thermodynamics,  tor  an  inviscid  and 
adiabatic  fluid  under  the  assumptions  of  irrotationality  as 
well  as  that  of  potential  flow,  the  governing  equation  for 
the  unsteady  aerodynamic  theory  is  given  by  a  form  of  tne 
wave  equation: 


1  3  2  $  3 ( V 2 )  V2 

7 2  $  -  -  [  -  +  -  +  V  •  V  (  -  )  ]  =0 

a2  3  1 2  3  1  2 


[2.  1 ) 


where  $  is  the  velocity  potential,  V  is  the  flow  velocity 


vector,  and  a  is  the  speed  of  sound. 


Assuming  harmonic 
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motion  and  using  linear  superposition  of  singularities  on 
the  wing,  the  solutions  to  this  equation  can  be  written  in 
integral  form  as  a  relationship  between  the  velocity  normal 
to  the  wing  and  the  pressure  on  the  wing. 

Two  basic  classes  of  lifting  surface  theories  now 
widely  used  are  the  kernel  function  method  in  its 
traditional  form  and  the  doublet  lattice  method  wmch 
represents  a  finite  element  technique. 

In  the  kernel  function  method,  assumed  pressure 
functions  are  required  as  input  to  these  solutions  and  the 
unknowns  are  the  coefficients  associated  with  these  pressure 
functions. 

The  finite  element  methods  io  not  require  assumed 
pressure  functions  as  input  to  obtain  a  solution.  This 
advantage,  however,  is  offset  by  the  tact  that  the  cost  per 
solution  is  much  higher,  as  the  finite  element  method 
requires  a  very  fine  mesh  and  thus  a  large  number  of 
elements  to  converge  to  an  acceptable  solution. 

On  the  other  hand,  if  special  weighting  functions  are 
used  in  the  assumed  pressure  distributions  with  implicit 
discontinuous  characteristics,  the  kernel  function  method 
can  show  a  much  faster  rate  of  solution  convergence  than  the 
finite  element  method.  It  is,  thus,  obvious  that  the  kernel 


function  method  can  have  an  economic  advantage  over  finite 


element  methods. 


■'.'HA 
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To  investigate  an  F-5  wing  (Figure  1)  with  a  tip 
mounted  launcher  and  AIH-7S  missile  (Figure  2)  in  unsteady 
flow,  the  studies  herein  were  conducted  by  using  a  kernel 
function  code  (ANKF)  for  subsonic  and  supersonic  flow  and  a 
doublet  lattice  code  (H7WC)  tor  only  subsonic  flow. 

2.1  DOtJBLET  LATTICE  SET HOD 
2-1.1  Introduction 

To  account  for  the  interference  between  a  body  (sucn  as 
a  fuselage)  and  a  wing  and  between  a  wing  (such  as  a 


stabilizer)  and  a  wing,  lifting  surface  finite  elements 
which  are  not  necessarily  planar  are  applied  to  both  the 
body  surfaces  and  lifting  surfaces  in  this  method.  Here, 
any  body  is  regarded  as  an  annular  wing  as  was  used  by 
Woodward  (13)  in  his  study  ot  the  steady  case.  This  method 
calculates  lift  distributions  on  surfaces  in  steady  and 
oscillatory  motion,  even  in  the  case  of  nonplanar  surface 
configuration,  at  only  subsonic  speeds. 

The  doublet  lattice  method  is  simple,  accurate,  and 
very  versatile.  An  overview  of  the  developments  of  the 
steady  and  unsteady  doublet  lattice  methods  is  given  in  the 
paper  by  Giesing,  Kalman,  and  Rodden  (20). 
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UNITS  IN  INCHES 


2.1.2 


Theoretical  Development 

The  forces  exerted  on  a  wing  depend  on  the  velocity 
components  of  the  air  stream.  These  are  the  velocity  or  tne 
free  stream  U m  which  is  opposite  to  the  direction  of  motion, 
and  the  velocity  w  which  is  perpendicular  to  the  direction 
of  motion.  Here,  the  velocity  normal  to  an  oscillating 
surface,  W,  can  be  calculated  from  the  equation 


M  =  U  Re(we1a)t) 


(2.2) 


w  here 


w  ( x  ,y ,  z  )  =  —  K(x-5,y-n.z-c.u>,M)AC_d£da 

8v  p 

L.S. 


(2.3) 


where  £  is  the  streamwise  coordinate,  a  is  the  tangential 
spanwise  coordinate  (Figure  3),  and  w  is  tne  frequency  of 
the  oscillation.  K  ()  is  the  kernel,  1  is  the  dach  nuaaer, 
and  AC p  is  the  lifting  pressure  coefficient  which  can  be 
expressed  in  terms  of  lifting  pressure  and  dynamic  pressure 


as  follows: 


e 
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In  ♦■ho  integration,  the  functions  are  integrated  over  all  oi 
the  lifting  surfaces. 

Since  the  doublet  lattice  method  applies  the  lifting 
surface  finite  elements  over  all  the  nody  surfaces  and  tae 
wing  surfaces,  the  above  integration  can  be  written  as 
follows : 


w(x,y,z)  =  — EAC 

8tt  ps 


j  K(x-£,y-n,z-s,a>,M)d£;da 

El ement 


(2.4) 


where  subscript  s  indicates  the  sending  elements  of  the 
surfaces. 

The  basic  concept  of  the  doublet  lattice  method,  which 
determines  the  influence  of  an  oscillating  lifting  surface 
element  at  a  point,  assumes  that  the  lifting  pressure  could 
be  concentrated  along  a  line.  This  line  is  located  at  cne 
1/4  chord  line  of  the  sending  element  (Figure  4).  Here,  the 
steady  flow  effect  is  modeled  by  a  horseshoe  vortex  which  is 
located  in  each  element  with  its  bound  leg  along  this  1/4 
chord  iine.  Additionally,  the  unsteady  flow  effect  is 
modeled  by  a  uniform  line  of  acceleration  potential  doublets 
which  is  superimposed  onto  the  bound  vortex.  Thus,  both 
vortex  and  doublet  effects  are  represented  by  an  influence 
function  which  relates  the  pressure  amplitude  and  the 
norma lwash  amplitude  at  a  given  control  point.  Therefore, 
the  first  step  of  the  procedure  reguires  the  dividing  of  the 
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Figure  4:  Idealization  ot  Wing  Panel  into  3oxes 
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lifting  surfaces  into  snail  boxes  (elements)  which  are 
oriented  in  strips  parallel  to  the  free  stream.  Moreover, 
it  is  advantageous  to  align  any  surface  edges,  hinge  lines, 
and  told  lines  with  the  boundaries  of  the  elements.  Tnen, 
the  integral  equation  (2.4)  can  be  rewritten  as  follows: 

AC  r 

w(x  ,y  ,z)  =  ZAC  --  K(x-c  i  ,  ,y-n » z-c  »u) ,M )do  (2.5) 

PS8tt  /4 

El ements 

In  this  eguation,  the  normalwash  boundary  condition  w(x,y,z) 
is  known,  but  the  lifting  pressure  coefficient  ACp  over 
each  element  is  unknown. 

James  (21)  used  a  control  point  (or  receiving  point)  at 
which  the  normalwash  boundary  condition  is  satisfied  on  each 
element.  He  located  this  point  at  the  center  in  spanwise 
direction  on  the  J/4  chord  line  of  each  element.  For  a 
two-dimensional  case,  he  showed  that  this  location  gave 
highly  accurate  results  and  was  considered  to  be  optimum  for 
the  doublet  lattice  method. 

Therefore,  in  order  to  calculate  the  unknown  values  of 
ACp  ,  we  can  write  equation  (2-5)  in  matrix  form: 

{  w  }  »  [  D  ]  {  ACp  }  (2.6) 


where  the  elements  of  the  normalwash  factor  matrix  [  D  ]  are 
given  as 


D  =  K(x-c1^,y-n,z-c,cj,M)da  (2.7) 

El ement$ 

In  general,  we  solve  this  matrix  equation  by  inverting  the 
matrix  [  D  ]  as  follows: 


{  ACp  }  =  [  D  ]_1  {  w  } 


(2.8) 


However,  when  the  matrix  [ 0  ]  is  large,  it  is  not  easy  and 
not  efficient  to  calculate  its  inverse.  Therefore,  tne 
present  method  solves  the  linear  simultaneous  equations 
(2.6)  by  direct  Gaussian  Triangular ization  and  uack 
solution.  Due  to  this  fact,  the  present  doublet  lattice 
method  has  the  disadvantage  that  additional  problems  can 
only  be  solved  by  repeating  the  whole  procedure  of 
triangulariza tion  and  back  solution. 

2.2  KERNEL  FUNCTION  METHOD 
2.2.1  Introduction 

While  the  doublet  lattice  method  can  only  calculate  the 
subsonic  flow  case,  the  kernel  function  method  can  be 


applied  not  only  to  the  uniform  and  non-uniform  subsonic 


flow  c< 153,  but  also  to  the  transonic,  and  supersonic  fiow 
cases.  Moreover,  we  can  apply  local  Macn  number 
distributions  over  the  lifting  surfaces  in  this  method,  ana 
include  the  dynamic  effects  of  the  components  of  the  lifting 
surfaces  and  the  external  stores. 

This  method  uses  assumed  pressure  functions  to  solve 
the  linear  lifting  surface  problem  with  interference  in 
steady  or  oscillating  flow.  Sines,  in  the  assumed  pressure 
functions,  a  special  weighting  function  is  used  to  imply  tne 
discontinuous  characteristics  of  the  pressure  distributions, 
the  solution  converges  faster  than  is  possible  with  the 
conventional  pressure  functions  composed  of  smooth 
functions.  This,  however,  also  presents  the  greatest 
difficulty  in  the  use  of  this  method  since  proper 
assumptions  have  to  be  made  for  selecting  the  pressure 
functions.  Details  of  the  kernel  function  method  are 
presented  by  Cunningham  in  deference  id  . 

2.2.2  Theoretical  Development 

As  in  the  doublet  lattice  method,  the  forces 
experienced  by  the  surfaces  are  related  to  the  downwash 
which  is  the  velocity  normal  to  the  lifting  surface.  The 
relationship  between  downwash  and  pressure  can  be  expressed 
as  follows: 


w ( 1 ,m,n )  = 


-  E  [  AP (a  ,  8  ,  y  )  K ( 1 -a ,m-8  »  n-Y  »  k ,M ) d! 
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whore  the  downwash  w(l,ra,n)  is  the  normal  velocity  at 
control  point  (l,m,n),  and  the  function  AP'a,B,y)  is  the 
normal  pressure  iistributior.  on  the  lilting  surface  at  a 
load  (or  integration)  point.  K  ()  is  the  kernel  function,  k 
is  the  reduced  frequency,  and  M  is  the  Mach  number. 

In  equation  (2.9),  while  the  downwash  w(l,m,n)  is  a 
known  quantity  and  is  given  by  the  boundary  conditions  on 
the  respective  lilting  surfaces  as  depending  on  the  mode 
shape  and  the  frequency  of  oscillation,  A?(a,B,y)  is 
unknown.  Therefore,  to  solve  this  equation  by  the  present 
method,  the  unknown  pressure  functions  are  assumed  to  be  a 
set  of  chordwise  and  spanwise  polynomials  with  weighting 
functions.  The  form  of  the  assumed  pressure  functions  is  as 
follows: 


AP  (a  ,  0  )  =  qPw(a,6)  [  90  (  §  )^0<« )  +  +  •••  ]  (2.10) 

where  the  coordinates  a  and  jjj  are  defined  in  the  transformed 
plane  of  the  surface  (Figure  5),  and  g  is  the  freestream 
dynamic  pressure.  Here,  the  chordwise  functions,  f  n  ( ci )  ,  are 

y;> -  uo<“>  - 1 

f-|(a)  =  U-|(a)  +  =  2a  +  1 

fn(«)  =  un(5)  +  Vi(“J 


(2.  11) 
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and  the  spanwiso  functions,  0^(3),  are  expressed  as 


q  (q)  =  [  a  U  ( 3 ) 
ym'  '  L  mo  o'  ' 


+  amlUl(i) 


(2. 12) 


where  the  coefficients  "a"s  are  the  unknowns  and  I  (  )  are 

n 

the  Tschetychev  polynomials  of  the  second  kind: 

U0(§)  »  1 
11,(5)  -  26 

.  (2.13) 

un(5)  =  2eun_n (e)  -  un_2( b ) 

In  equation  (2.10),  the  function  P  ( ci,  g)  is  the 

w 

weighting  function  which  depends  on  the  pressure  types  in 
both  chordwise  and  spanwise  directions  as  shown  in  Figure  t> , 
and  on  the  speed  of  the  flow  at  the  downwash  point. 

If  the  downwash  point  is  subsonic,  a  subsonic  weighting 
function  is  used.  Similarly,  if  the  downwash  point  is 
supersonic,  a  supersonic  weighting  function  may  be  used  to 
calculate  the  downwash  and  the  interference  effect.  In  this 
report,  the  function  P  (a,S)  was  assumed  as  follows: 

W 


4PU  2 

pw(5.5)  =  — t~~ 


3 


S  h ( a ) 1 (8) 


(2. 14) 


22 


and  was  used  regardless  of  Mach  number  for  subsonic, 
transonic,  and  supersonic  solutions.  Additional  discussions 
about  a  supersonic  weighting  function  are  given  in  Reference 
18  . 
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Chapter  III 
NOHEBICAL  HOBBLING 


In  order  to  model  the  wing  with  tip  mounted  launcher 
and  missile,  it  was  assumed  that  the  wing  and  launcher  be 
flat  plates.  The  dimensions  of  the  wing  were  the  same  as  in 
the  NLR  experiments  where  a  slightly  modified  haltmodel  of 
the  outer  part  of  the  P-5  wing  [scale  1:4.5)  was  used.  The 
missile  was  simplified  without  fins  as 

(1)  a  flat  plate  for  accuracy  evaluation 

(2)  an  endplate 

(3)  various  flat  plate  -  endplate  combinations 
to  account  for  cross  flow  and  pressure  drop 

(4)  an  octagonal  cylinder  open  at  both  ends 

(5)  an  axisymmetric  slender  body 

[for  the  doublet  lattice  method  only) 

For  case  3,  various  cross  sections  were  considered 
[Figure  7),  and  tor  case  4,  a  cylinder  was  idealized  as 
having  an  octagonal  shape  with  eight  lirting  surfaces.  It 
is  obvious  that  case  4  needs  a  large  number  of  elements  thus 
increased  the  computer  cost. 

For  both,  doublet  lattice  and  kernel  function 
computational  methods,  symmetric  flow  with  respect  to  the 
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Figure  7: 


Various  Tip  Store  Cross  Sections 
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x- ixis  of  the  wing  was  assumed  and  the  dach  number  was 
assumed  to  be  uniformly  distributed  over  the  wing  surface. 

To  calculate  the  vibration  modes  of  the  wing  and  store, 
the  following  polynomial  expression  was  taken  from  Reference 
9: 


wing  :  m= 0,1,2  n=0,l 

w(x,y)  =  I  amnxV 

store  :  m=0  n=0,l 

This  expression  assumes  linear  displacements  in 
chord  wise  direction  and  parabolic  displacements  in  span wise 

direction.  The  coefficients  a  are  shown  in  Table  1  ror 

tnn 

several  dach  numbers.  In  the  experiments,  the  wing  model 
was  oscillated  in  pitch  about  a  50?  root  chord  axis  at 
varying  frequencies. 

A  sample  mode  shape  for  the  clean  wing  as  used  for  the 
computation  of  the  normalwash  boundary  conditions  in  the 
unsteady  flow  case  at  d=0.95,  F=20Hz,  ind a =0  is  exhibited 

in  Figure  8  .  The  i isplacemeu ts  are  magnified  to  better 
show  the  deflection  shape. 

In  this  paper,  a  dach  number  range  between  0.6  and  1.35 
is  covered  while  the  frequency  or  oscillation  of  the  wing  is 
kept  at  20Hz.  generally,  when  the  frequency  is  increased. 


the  imaginary  component  of  the  lift  increases  but  the  real 
Dirt  of  the  lift  decreases  very  slightly  at  constant  dach 


number.  When  the  Mach  number  is  increased  without  varying 
the  frequency,  the  real  component  of  the  lift  is  increased 
while  the  imaginary  component  is  decreased. 

It  was  assumed  that  the  lower  surface  and  upper  surrace 
of  the  wing  have  the  same  aerodynamic  characteristics  for 
simplicity.  The  geometric  characteristics  of  the  wing  are 
as  follows: 


Reference 

Area 

= 

0.  26045 

m2 

f 

403.7 

in2 

Reference 

Chord 

= 

0.4183 

a 

( 

16.4685 

in 

Reference 

Seoispan 

— 

0.6226 

m 

( 

24.5118 

in 

Wing  Section  %  Span 

1  18.1 

2  35.2 

3  51.2 

4  64.1 

5  72.1 

6  81.7 

7  87.5 

8  97.7 

Here  the  wing  sections  were  chosen  the  same  as  those  of 


the  MLR  experiments  in  order  to  compare  computational  and 
experimental  results. 
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3.1  DOOBLBT  LATTICE  METHOD 

The  model  used  for  the  doublet  lattice  method  is  shown 
in  Figure  9  .  The  wing  is  formed  by  11  strips  which  are 
subdivided  into  10  chordwise  elements  or  boxes.  Tne 
launcher  is  idealized  by  2  strips  which  have  lb  boxes  or 
elements  in  chordwise  direction  each.  Since  our  main 
concern  is  the  wing  tip  section,  more  elements  have  been 
located  in  that  area.  The  tip-missile  model  is  shown  in 
Figure  10  .  The  flat  panel  which  represents  the  missile 
model  consists  of  2  strips  with  a  total  of  19  boxes  each. 
The  width  of  the  plate  for  the  missile  model  is  equal  to  tne 
diameter  of  the  actual  missile.  For  this  method  only,  a 
cylindrical  axisymaetric  slender  body  model  was  also  used. 

It  should  be  mentioned  here  again  that  this  method  can 
only  be  applied  to  subsonic  flow. 

3-2  KERNEL  FUNCTION  METHOD 

The  model  used  tor  the  kernel  function  method  is 
similar  to  the  doublet  lattice  model.  Instead  of  boxes,  it 
has  structural  points  at  which  mode  deflections  are  given. 
There  are  83  structural  points  on  the  wing,  19  structural 
points  on  the  launcher,  and  11  structural  points  on  the 
interface  of  wing  and  launcher  (Figure  11). 

The  modeling  of  the  tip-missile  is  almost  the  same  as 
for  the  doublet  lattice  case.  There  are  21  structural 
points  on  the  flat-panel  modeled  missile,  and  18  interface 


s 
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structural  points  between  launcher  and  tip-missile  (Figure 
12).  For  this  method,  supersonic  cases  can  be  analyzed 
additionally  but  the  slender  body  model  cannot  be  used. 

The  spanwise  loading  ot  the  wing  was  assumed  to  be  ot 
type  2  in  Figure  6  and  the  surface  to  be  of  type  1.  For  the 
launcher,  a  spanwise  loading  of  type  2  and  the  surface  type 
3  were  used.  Various  different  surface  and  loading  types 
were  used  for  the  missile  depending  on  the  shape  of  the 
model.  The  chordwise  loadings  of  wing,  launcher,  and 
missile  were  assumed  to  be  of  type  1. 


SLENDER  BODY 
MODEL 


OCTAGONAL  CYLINDER 
OPEN  AT  BOTH  ENDS 


Figure  10: 


geometrical  .lodel  of  dissile  'DL'^) 


LAUNCHER  STRIP 


FLAT  PLATE 


LAUNCHER  L.E. 


WING  L.E.  WING  T.E.  LAUNCHER  T.E 


TOP  VIEW 


Figure  12:  Geonetrical  :iodei  of  *.issile  'K?X) 
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BESULTS 


In  all  cases,  the  wring  was  oscillated  about  a  near  zero 
mean  angle  of  incidence.  Thus,  the  unsteady  results  are 
based  on  the  total  unsteady  pressure  difference  between 
upper  and  lower  surfaces. 

Only  curves  for  the  lift  coefficient  distributions  over 
the  wing  are  presented  since  the  moment  coefficient 
distributions  depend  on  the  location  of  the  elastic  axis  of 
the  wing  and  can  be  computed  by  direct  integration  of  tne 
pressures  over  the  wing.  Lift  and  pressure  coefficients  for 
the  various  store  models  were  obtained  but  are  not  given 
here  since  they  are  very  small  and  vary  widely  between  the 
different  models.  Also,  no  steady  results  are  presented, 
since  the  assumption  of  the  wing  as  a  flat  plate  without 
camber  and  thickness  does  not  yield  agreeable  results. 


4.1  CLB1M  BING  CASE 

Before  analyzing  the  characteristics  ot  the  various 
tip-store  models,  it  is  necessary  to  verify  noth  computer 
codes  by  computing  unsteady  clean  wing  results  and  comparing 
these  to  existing  results  such  as  those  of  Parker  and  the 
experimental  results  of  NLR. 

The  clean  wing  was  analyzed  in  subsonic  flow,  Jl=0.6, 
and  at  a  frequency  of  oscillation  of  F=20Rz.  As  evident  in 
Figure  1J  and  Table  2,  the  unsteady  lift  coefficient 
distributions  from  the  doublet  lattice  method  (DLM)  and  from 
Parker  exhibit  good  agreement,  especially  at  the  root  and 


tip  sections  ot  the  wing. 


It  is  believed  that  tne 


discrepancy  between  the  present  DLM  and  Parker  at  60%-90a  of 
the  span  of  the  wing  might  have  resulted  from  different 
vibration  mode  coefficients. 

The  results  from  the  kernel  function  method  show  a 


slightly  different  shape  than  the  two  others. 


Here,  the 


lift  coefficient  distribution  presents  higher  values  than 
the  other  two  methods  from  the  root  section  to  mid-span. 
From  mid-span  to  the  tip-section,  the  lift  coefficient  drops 
off  faster  than  for  the  other  methods,  but  it  is  still  close 
to  both,  the  present  DLM's  and  Parker's  results. 

All  three  sets  of  results  show  values  higher 
than  the  experimental  results  from  NLR.  This  migat.  be  due 
to  the  fact  that  no  fuselage- w ing  interaction  is  considered 
nor  is  the  wing-tip  upwash  condition  being  accounted  for. 
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TABLE  1 


■ 

Coefficients 

"a_  ••  for 
mn 

V  arious 

Vibration 

loies 

M=0.  6 

M=0.  7 

.'1=0.3 

3=0.9 

1 

k  =  0. 200 

k=0. 173 

k=0. 163 

k  =  0. 133 

aOO 

-0. 336 

-0.346 

-0. 351 

-0.349 

,  • 

aO  1 

0.  997 

1.020 

1. 022 

1.016 

v*. 

Wing  alO 

0.  045 

0.127 

0.  162 

0.  082 

all 

0.032 

-0.260 

-0.296 

-0.209 

m 

a  20 

-0.141 

-0.  40  0 

-0. 531 

-0. 442 

a  2 1 

-0. 001 

0.596 

0.677 

0.424 

— 

Tipstore 

aOO 

-0. 389 

-0.450 

-0.459 

-0.513 

aO  1 

1.  022 

1.099 

1. 075 

1.099 

TABLE  2 

Unsteady  Lift  Coefficient  Distributions  of  Clean  Wing 


Wing  Span 
(*) 

Parker 

d.l.  a 

Error 

:*> 

ANK? 

Error 

f  w  \ 

[  o) 

9.  75 

1.04 

1.0  4 

0.0  0 

1.03 

3.35 

28.  89 

1.16 

1.  1  7 

0.  36 

1.20 

3.45 

46.  91 

1. 27 

1.23 

0.  79 

1.31 

3.  1  5 

63.  13 

1.36 

1.39 

1.4/ 

1.3  7 

0.  74 

76. 92 

1.  38 

1.4  1 

2.  1  7 

1.38 

0.00 

9  7. 76 

1.30 

1.34 

3.03 

1.24 

-4.52 

95.  23 

1.06 

1.09 

2.33 

0.8  3 

-21. 70 

99. 04 

0.67 

0.69 

2. 93 

0.4  9 

-26.37 

EXPERIMENT 


Since  the  DLM  shows  only  about  3%  difference  with 


I 

a 


Parker's  results,  this  code  seems  promising, 
though  the  difference  between  the  results  from 
function  code  and  from  Parker  is  somewhat  large 
code  shows  still  good  enough  comparison  to  be 
analysis  since  its  distribution  is  clos 
experimental  results  near  the  wing  tip. 


And  ever* 
the  kernel 
r,  the  ANKf 
used  for  the 
er  to  the 
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4.2  DOOBLBT  LATTICE  METHOD 

The  22  models  mentioned  previously  were  investigated  i 
the  subsonic  and  transonic  region  for  Mach  numbers  fro 
M  =  0.6  -  0.9  by  use  of  the  doublet  lattice  method  cod 
(H7WC) .  Plots  of  the  unsteady  lift  coefficien 
distributions  are  shown  in  Figures  14  through  17  for  th 
wing  with  the  missile  modeled  by  different  combinations  o 
plates,  by  an  octagonal  cylinder,  and  by  a  cylindrica 
slender  body. 

As  shown  in  the  figures,  the  results  of  the  DLM  presen 
higher  distributions  than  the  experimental  results 
Compared  to  the  clean  wing,  the  missile  increases  the  load 
over  the  entire  wing  with  a  greater  effect  in  the  ti 
region.  As  can  be  seen,  the  lift  coefficient  distribution 
of  the  endplate  models  show  better  comparison  with  th 
experimental  result  in  the  tip  region.  Model  2,  where  th 
middle  of  the  endplate  is  positioned  next  to  the  launcher 
reveals  very  good  agreement  with  model  iy,  where  th 
endplate  is  located  above  the  outer  launcher  edge. 


L 

m 

e 

t 

e 

f 

1 

t 


t 

£ 

e 

e 

t 

e 


\ 


i 


L 


Figure  15: 


Onsteady  Lift  Coefficient  Distrioutions  ' 
k  =  0.  173,  DL:i,  22  Models) 
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Another  endplate  model,  model  20,  where  the  endplate 
lies  underneath  the  outer  launcher  edge,  shows  a  higher  lift 
distribution  in  the  tip  region.  Model  1,  m  which  the 
missile  was  modeled  by  a  flat  plate,  and  model  1 8,  where  the 
missile  was  modeled  by  three  plates  at  120  degrees  to  each 
other,  show  distributions  similar  to  each  other.  As  stated 
earlier,  the  flat  plate  model  is  widely  used  at  present  tc 
model  external  bodies  or  stores. 

The  cylindrical  slender  body,  model  22,  shows  the 
lowest  loads  in  the  tip  region.  The  other  models  such  as 
the  combinations  of  flat  plates  and  endplates,  the  octagonal 
cylinder,  and  the  cruciform  model  presented  increased  loads 
in  the  tip  region.  The  loads  distribution  tor  these  models 
can  hardly  be  differentiated.  Therefore,  to  compare  the 
models  in  more  detail,  8  representative  models  were  picxed 
out  of  the  previous  22. 

1.  Flat  plate  'Model  1) 

2-4.  Endplate  -  J  different  types  (Models  2,  19,  and  20) 

i>.  Combination  of  flat  plate  and  endplate  (Model  4) 

6.  Cruciform  (Model  14) 

7.  Octagonal  cylinder  (Model  21) 

8.  Slender  body  cylinder  (Model  22) 

The  unsteady  lift  coefficient  distributions  of  these  eight 
models  together  with  experimental  results  at  M=0.6  are 
presented  in  Tables  J  and  4  and  in  Figure  18  . 


MODEL  4 
MODEL  14 
MODEL  19 
MODEL  20 
MODEL  2 
MODEL  22 
EXPERIMENT 
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Unsteady  Lift  Coefficient  Dintr ibu tions  ':i=0.6 
k  =0.  20  0  ,  DLM,  3  Kod  els) 


The  characteristic  bump  present  in  the  experimental 
litt  distribution  at  about  60*  semispan  is  barely  visible  in 
the  computational  results  and  is  shitted  to  about  70/t 
semispan.  Overall,  with  the  exclusion  of  this  local 
increase  in  lift  coefficient,  the  best  computational  results 
were  those  tor  the  cylindrical  slender  body  model  and  were 
about  10*  larger  than  the  experimental  data  across  the  winy 
and  about  l5%-20%  larger  in  the  tip  region. 

When  the  flow  speed  was  increased  to  M=0.7  and  the 
reduced  frequency  thus  lowered  to  k=0.1bd,  lift  was 
increased  over  the  entire  wing  (Figure  19,  Tables  b  and  b) . 
Also,  the  eight  models  show  similar  distributions  as  tor 
f1=0 . 6 .  The  difterences  to  the  experimental  results, 
however,  are  increased. 

As  the  flow  speed  increases  to  11=0. b  and  11  =  0.9,  tae 
litt  coefficients  ot  all  models  are  increased.  Results  are 
presented  in  Figures  20  and  21  and  in  Tables  7  through  10  . 
In  these  linearized  transonic  cases,  the  cylindrical  slender 
body  model  shows  results  closest  to  the  experimental  data 
together  with  the  endplate  models  while  the  cruciform  and 
octagonal  cylinder  models  display  increased  lift  in  the  tip 
region.  The  model  with  the  flat  plate  -  endplate 
combination  presents  the  highest  lift  distribution  over  the 
entire  wing  for  every  case.  Since  the  linear  code  cannot 
predict  the  nonlinear  transonic  flow,  results  do  not  snow 
good  comparison  to  experimental  data. 


TABLE  5 


Unsteady  Litt  Coefficient  Distriuutioas  (Real,  M=0.7, 
k=0. 173,  DIM,  8  Models) 


Wing 

Section 

1 

2 

3 

4 

5 

6 

7 

8 

Model 

1 

1.21 

1.3b 

1.51 

1.65 

1.74 

1.8  1 

1. 84 

1.79 

Model 

2 

1.  20 

1.34 

1.-9 

1.62 

1.70 

1.  75 

1.77 

1.  o7 

Model 

4 

1.22 

1. 37 

1.53 

1.67 

1.77 

1.86 

1.91 

1.  92 

Model 

1  4 

1.22 

1.36 

1.52 

1.66 

1.75 

1.83 

1.87 

1. 85 

Model 

19 

1. 20 

1.34 

1.48 

1.61 

1.70 

1.75 

1.  lb 

1.66 

Model 

20 

1.21 

1.35 

1.50 

1.63 

1.72 

1.78 

1.81 

1.  7  o 

Model 

2  1 

1.21 

1.36 

1.51 

1.65 

1.74 

1.82 

1.  86 

1.83 

Model 

22 

1.  19 

1.33 

1.47 

1.59 

1.67 

1.70 

1.70 

1.50 

SXP 

1.00 

1.11 

1.23 

1.29 

1.24 

1.3  1 

1.  29 

1.15 

TABLE  6 

Unsteady  Litt  Coefficient  Distributions  (Imag,  M=0.7, 
*=0.173,  DLM,  8  Models) 


Wing 

Section 

1 

2 

3 

4 

5 

6 

7 

8 

Model 

1 

0.35 

0.  35 

0.34 

0.34 

0.33 

0.32 

0.31 

0.28 

Model 

2 

0.36 

0.  35 

0.35 

0.  34 

0.34 

0.32 

0.31 

0.27 

Model 

4 

0.35 

0.  35 

0.  34 

0.33 

0.33 

0.3  1 

0.31 

0.29 

Model 

14 

0.  35 

0.  35 

0.  34 

0.34 

0.33 

0.32 

0.31 

0.29 

Model 

19 

0.  36 

0.  35 

0.  35 

0.  34 

0.33 

0.  3  1 

0.  30 

0.26 

Model 

20 

0.  36 

0.35 

0.  35 

0.34 

0.34 

0.32 

0.  31 

0.  29 

Model 

21 

0.  35 

0.  35 

0.34 

0.33 

0.33 

0.3  1 

0.30 

0.27 

Model 

22 

0.36 

0.  36 

0.35 

0.34 

0.34 

0.  32 

0.  30 

0.25 

EXP 

0.24 

0.  19 

0.  1b 

0.12 

0.12 

0.10 

0.  09 

0.09 

TABLE  7 

Unsteady  Litt  Coefficient  Distributions  (Real,  M  =  0.8 
k=0. IbJ,  DLM,  8  Models) 


Ming 

Section 


Model 

Model 

Model 

Model 

Model 

Model 

Model 

Model 

EXP 


b 

b 

7 

8 

.89 

1.9b 

2.  00 

1.94 

1 

.84 

1.90 

1.92 

1.81 

'•  j 

.92 

2.02 

2.07 

2.07 

.90 

1.99 

2.  03 

1.99 

.84 

1.90 

1.91 

1.  80 

.8b 

1.93 

1.9b 

1.90 

.89 

1.98 

2.  02 

1.99 

** 

.80 

1.84 

1.84 

1.  62 

.29 

i.b  7 

1.44 

1.09 

j 

TABLE  8 


Unsteady  Lift  Coefficient  Distributions  (Imag,  3=0.8 
k=0.1b3,  DLM,  8  Models) 


Model  1 
Model  2 
Model  u 
Model  14 
Model  19 
Model  20 
Model  21 
Model  22 


2 

3 

.  28 

0.27 

.29 

0.27 

.27 

0.  26 

.  28 

0.  27 

.29 

0.27 

.28 

0.27 

.  28 

0.  26 

.29 

0.  28 

.  23 

0.16 

.2b 

0.2b 

.26 

0.26 

.2b 

0.24 

,2b 

0.2b 

.26 

0.2b 

.26 

0.26 

.2b 

0.24 

.27 

0.26 

TABLE  9 


l 

i 

Unsteady 

Lift  Coefficient  Distributions 

(Real 

,  M  =  0 . 

9, 

k=0. 1 38 

,  DIM, 

8  Models) 

L 

Wing 

1 

2 

3 

4 

5 

6 

7 

8 

* 

Section 

Model 

1 

1.44 

1.61 

1.77 

1.92 

2.01 

2.06 

2.  09 

2.01 

Model 

2 

1.43 

1.58 

1.  74 

1.88 

1.96 

1.99 

2.  00 

1.86 

Model 

4 

1. 4b 

1.62 

1.80 

1.95 

2.0  5 

2.  13 

2.  17 

2.15 

AC 

Model 

1  4 

1. 45 

1.61 

1.78 

1.93 

2.03 

2.09 

2.12 

2.07 

Model 

19 

1.42 

1.58 

1.74 

1.88 

1.96 

1.99 

1.99 

1.86 

Model 

20 

1. 43 

1.  59 

1.76 

1.90 

1.98 

2.03 

2.0  4 

1.  9b 

Model 

21 

1. 44 

1.61 

1.  78 

1.92 

2.02 

2.08 

2.  1 1 

.  2.06 

Model 

2  2 

1.41 

1.57 

1.72 

1.85 

1.91 

1.92 

1.91 

1.67 

EXP 

0.87 

1.38 

1.  87 

1.84 

1.45 

1.61 

1.77 

1.66 

u 

11 

TABLE 

10 

Unsteady 

Lift  Coefficient  Distributions 

(Imag 

• 

o 

It 

*3 

% 

9, 

: 

k=0. 1 38 

,  DLfi , 

8  Models) 

J 

Wing 

1 

2 

3 

4 

5 

6 

7 

a 

' 

Section 

1 

Model 

1 

0.20 

0.18 

0.  16 

0.13 

0.  12 

0.09 

0.  07 

0.05 

Model 

2 

0.21 

0.  20 

0.17 

0.15 

0.13 

0.10 

0.09 

0.  0  t> 

8 

Model 

4 

0.  19 

0.  17 

0.  15 

0.  12 

0.  1  0 

0.07 

0.  06 

0.  04 

\ 

Model 

14 

0.19 

0.18 

0.  15 

0.13 

0.  1  1 

0.08 

0.  07 

0.05 

\  - 
»*.  • 

Model 

19 

0.21 

0.  19 

0.17 

0.15 

0.13 

0.1  0 

0.  08 

0.05 

* 

Model 

20 

0.21 

0.19 

0.  17 

0.  14 

0.  13 

0.  1  0 

0.  09 

0.07 

Model 

21 

0.  19 

0.  18 

0.  15 

0.  12 

0.  1 1 

0.07 

0.  06 

0.03 

Model 

22 

0.22 

0.21 

0.  1  8 

0.  16 

0.1  5 

0.1  1 

0.  10 

0.07 

EXP 

0.28 

0.13  - 

0.11- 

0.25 

-0.05 

-0.07 

-0.  27 

-0.  19 

-  MODEL  1 

-  MODEL  2 

-  MODEL  4 

-  MODEL  14 

- MODEL  19 

-  MODEL  20 

-  MODEL  21 

-  MODEL  22 

o  o  a  EXPERIMENT 


. . . . . . . I . I . . . I . I 

0  0.1  0.2  0.3  0.«  0.5  0.5  0.7  0.5  0.9  1.0 

y/s 


e  2  1  : 


Unsteady  Lift  Coefficient  Distributions  ;d=U.9 
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In  optimizing  the  missile  models,  the  effect  of  cost 
was  an  important  factor.  Thus,  the  amount  of  computer  time 
and  storage  for  these  eight  models  on  an  I3M  J081D  mainframe 
computer  is  listed  in  Tables  11  through  14  .  As  can  be  seen 
in  the  tables,  the  octagonal  cylinder  model  required  the 
largest  amount  of  CPU  time  and  storage  followed  by  the 
combination  model  and  the  cruciform  model.  The  cost  for  the 
flat  plate,  endplate,  and  slender  body  models  are  less  than 
those  for  the  others.  In  addition,  these  models  snow 
considerably  better  agreement  with  the  experimental  lift 
distributions. 

4.3  KERNEL  FUNCTION  METHOD 

In  order  to  verify  the  advantage  of  the  end-plate 
models,  fe  models  were  then  tested  by  the  kernel  function 
code.  Here,  the  octagonal  cylinder  model  was  not  used  due 
to  the  highest  cost  and  the  slender  body  cylinder  model  was 
not  tested  since  this  type  of  geometry  cannot  oe 
accommodated  by  the  ANKF  code. 

Thus,  six  models  were  examined  by  use  of  this  code  as 
listed  below: 

1.  Flat  plate  (Model  i) 

2-4.  Endplate  -  3  different  types  (Models  2,  iy,  and  20) 
b.  Combination  of  flat  plate  and  end-plate  (Model  4) 
Cruciform  (Model  14) 
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TAELE  11 

Comparison  of  Effort  (M=0.6,  k=0.200#  DLH,  3  Models) 


CPU  Time 
(sec) 

Storage  I/O 
(K-oia) 

Disk  I/O 

;exc?) 

Model 

1 

83.36 

2338.2 

6119 

Model 

2 

88.64 

2388.8 

6119 

Model 

4 

121.85 

3902.3 

1  02  31 

Model 

14 

128.87 

3970.8 

10231 

Model 

19 

88.91 

2391.3 

6119 

Mode  1 

20 

91.83 

2412.8 

6119 

Model 

21 

190.85 

6158.9 

15710 

Model 

22 

66.04 

1839.4 

5081 

TABLE  12 

ComDarison  of  Effort  (M  =  0.7,  k=0.173,  DLM,  8  Models) 


CPU  Time 
Jsec) 

Storage  I/O 
(K-ain) 

Disk  I/O 
(EXCP) 

Model 

1 

83.99 

2345.3 

6129 

Model 

2 

89.34 

2  397.1 

6129 

Model 

4 

133.22 

4005.0 

10241 

Model 

14 

128. 60 

3969.0 

10241 

Model 

19 

88.  41 

2388.6 

6129 

Model 

20 

88.07 

2  385.5 

6129 

Model 

21 

191.93 

6170.6 

1  5720 

Model 

22 

66.03 

1 889.4 

50  81 

5b 


TABLE  13 


r 

Comparison 

ot 

Etfort  (M  =  0. 

8,  k=0.l53, 

DLM,  8  Models) 

— * 

» 

:=* 

CPU  Time 

Storage  I/O 

Disk  I/O 

L 

(sec) 

(K-min) 

(EXC?) 

a 

Model 

1 

83.07 

2337. b 

b  1  29 

Mode  1 

2 

87.70 

2  382.5 

6129 

Mode  1 

4 

119. 44 

3831.8 

10241 

-■> 

Model 

14 

132.22 

3995.1 

1  0241 

a 

Mo  de  1 

19 

88.10 

2386.3 

6129 

Model 

20 

87.98 

2385.0 

6129 

"T, 

■  ■ 

Model 

21 

191.73 

6 168.9 

1  5720 

i 

Model 

22 

65.86 

1 887.9 

5081 

.  -• 

5' 


The  six  models  were  evaluated  at  the  same  Mach  numbers 
and  reduced  frequencies  as  in  the  previous  section.  At 
M=0.6  and  k=0.200,  tive  of  the  models  show  lift 
distributions  somewhat  similar  to  those  for  the  DLM  except 
in  the  tip  region  where  the  drop  in  lift  coefficients  occurs 
closer  toward  the  wing  tip  (Figure  22,  Tables  1b  and  1b) . 
Two  of  the  endplate  models,  model  2  and  model  19,  show 
fairly  good  lift  distributions  compared  to  the  experimental 
results  while  the  other  endplate,  model  20,  shows  increased 
lift  in  the  tip  region  by  a  significant  amount.  In  tue 
previous  section,  the  flat  plate  model  showed  a  higher  lift 
distribution  m  the  tip  region  than  this  endplate  model 
while  here  this  situation  is  reversed.  The  combination 
model,  which  showed  the  highest  lift  in  the  tip  region  in 
the  DLM,  presents  almost  the  same  shape  for  the  KFM  results. 

It  is  of  interest  to  note  that  the  cruciform,  model  14, 
shows  highly  increased  lift  over  the  entire  wing  surface. 
It  seems  that  a  singularity  inherent  in  the  matrix  due  to 
the  symmetry  of  the  model  about  its  longitudinal  axis  may 
have  caused  this  excessive  result. 

When  the  Mach  number  was  increased  and  the  reduced 
frequency  thus  decreased,  lift  was  increased  as  was  seen  in 
the  previous  section.  The  lift  distributions  for  these 
cases  are  presented  in  Figures  23  through  25,  and  in  Tables 
17  through  22  .  With  the  increase  in  Mach  numoer,  the 
discrepancy  between  the  numerical  and  experimental  results 


TABLE  15 


5  fc 


Unsteady 

Lift 

Coefficient  Distributions 
k=0. 200 ,  KFM,  6  Models) 

(Real, 

M=0. 6, 

rrs 

Wing 

Span 

0-0 

<*) 

28.  2 

54.1 

75.6 

91.0 

99.0 

iSS 

** 

Model 

1 

1.06 

1.26 

1.44 

1.58 

1.62 

1.57 

Model 

2 

1 . 05 

1.24 

1.41 

1.52 

1.50 

1.38 

Model 

4 

1.07 

1.27 

1.45 

1.60 

1 .64 

1.  59 

.  • 

Model 

14 

1-24 

1. 50 

1.85 

2.39 

3.29 

4.2/ 

Model 

19 

1 . 06 

1.25 

1.42 

1. 52 

1.47 

1.30 

Model 

20 

1. 05 

1.25 

1.42 

1.57 

1 .63 

1 . 64 

y 

TABLE  16 

Unsteady  Lift  Coefficient  Distributions  (Iuag,  M=0.6, 
k=0.200,  KFM,  b  Models) 


Sing 

Span 

0.0 

(*) 

28.  2 

54.  1 

75.6 

91.0 

99.0 

Model 

1 

0.31 

0.31 

0.31 

0.  31 

0.31 

0.31 

Model 

2 

0.31 

0.31 

0.31 

0.31 

0.31 

0-31 

Model 

4 

0.31 

0.31 

0.31 

0.  31 

0.31 

0.31 

Model 

14 

0.28 

0.  28 

0.28 

0.  28 

0.28 

0.28 

Mode  1 

19 

0.  34 

0.34 

0.34 

0.34 

0.34 

0.34 

Mode  1 

20 

0.28 

0.  28 

0.28 

0.  28 

0.28 

0.  28 

60 


v 


:*> 

TABLE  17 


Unsteady 

Lift 

Coefficient  Distributions 
k=0.173,  KFM  ,  6  Models) 

'Real, 

M=0. 7, 

as 

Wing 

0.0 

28.2 

54.1 

75.6 

91.0 

99.0 

* 

Span 

<*) 

tt# 

Model 

1 

1.13 

1.34 

1.54 

1.71 

1 .77 

1.72 

MOdel 

2 

1 .12 

1.32 

1.51 

1.65 

1  .b3 

1.51 

Model 

4 

1.  14 

1.35 

1.55 

1.  73 

1.79 

1.74 

Model 

14 

1.37 

1 .67 

2.08 

2.  74 

3.79 

4.  90 

Model 

19 

1.  13 

1.34 

1.52 

1.65 

1  .b0 

1.43 

Model 

20 

1.11 

1. 32 

1.52 

1.69 

1  .78 

1.79 

TAELE  18 

Unsteady  Lift  Coefficient  Distributions  (Imag,  M=0.7, 
k= 0- 173,  KFM,  6  Models) 


Wing 

Span 

0.0 

CM 

28.2 

54.  1 

75.6 

91.0 

99.0 

Model 

1 

0.25 

0.25 

0.25 

0.  25 

0.25 

0.  25 

Model 

2 

0.25 

0.25 

0.25 

0.25 

0.25 

0.25 

Model 

4 

0.25 

0.25 

0.25 

0.  25 

0.25 

0.  25 

Model 

14 

0.20 

0.20 

0.20 

0.  20 

0.20 

0.  20 

Model 

19 

0.  30 

0.30 

0.30 

0.30 

0.30 

0.  30 

Mode  1 

20 

0.21 

0.21 

0.21 

0.  21 

0.21 

0.  21 

-  MODEL  1 

- MODEL  2 

-  MODEL  4 

-  MODEL  14 

-  MODEL  19 

-  MODEL  20 

□  DO  EXPERIMENT 
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TABLE  19 

Unsteady  Litt  Coefficient  Distributions  (Real,  M=0.8, 
k=0.153,  KFil,  6  Models) 


W  ing 
Span 

0.0 

(5) 

28.  2 

U“ 

45 

• 

75. b 

91.0 

99.0 

Model 

1 

1.23 

1.46 

1.68 

1.87 

1.92 

1.86 

Model 

2 

1.21 

1.44 

1.64 

1.79 

1.77 

1.62 

Model 

4 

1. 24 

1.48 

1.70 

1.89 

1 .94 

1.86 

Model 

14 

1.b3 

1.87 

2.33 

3.  04 

4.12 

5.  23 

Model 

19 

1.22 

1.46 

1.66 

1.  80 

1.76 

1. 57 

Model 

20 

1.21 

1.44 

1.65 

1.83 

1.91 

1. 89 

f.\;  ■ 

\  ; 

I—*-* 

m 

TABLE  20 

.1 

Unsteady 

Litt  Coefficient  Distributions 

(Imag,  M  =  0.8, 

4 

Ml 

k=0.153,  KFM,  6  Models) 

•m 3 

1 

V" 

L  - 

Wing 

0.0  28.2  54.1  75.6 

91.0  99.0 

'  *3 

n 

Span 

(*) 

TABLE  21 


Unsteady  Lift  Coefficient  Distributions  (Heal,  M=0.9, 
lc=0. 1  38 ,  KFM,  6  Models) 


King 

0.0 

28.  2 

64.  1 

75.6 

91.0 

99.0 

Span 

{%) 

Model 

1 

1.39 

1.67 

1.92 

2.  1 1 

2.  15 

2.  06 

Model 

2 

1.37 

1.63 

1.86 

2.01 

1  .97 

1.79 

Model 

4 

1.43 

1.71 

1.97 

2.  16 

2.1  8 

2.06 

Model 

14 

1.86 

2.27 

2.84 

3.  65 

4.82 

5.  98 

Mode  1 

19 

1. 40 

1.66 

1.90 

2.04 

1.97 

1.75 

Model 

20 

1.36 

1.62 

1.87 

2.05 

2.1  1 

2.  08 

TABLE  22 

Unsteady  Lift 

Coefficient  Distributions 
k=0.138,  KFM,  6  Models) 

(Iraag,  M=0.9, 

Wing 

0.0  28.2  54.1  75.6 

91.0  99.0 

Span 

(*) 

Model 

1 

0. 14 

0.14 

0.14 

0.  14 

0.14 

0.  14 

Model 

2 

0. 15 

0.16 

0.  15 

0.  16 

0.15 

0.  15 

Model 

4 

0.  14 

0.14 

0.  14 

0.14 

0.14 

0.  14 

Model 

14 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

Model 

19 

0.23 

0.23 

0.23 

0.  23 

0.23 

0.  23 

M  ^del 

20 

0.  06 

0.06 

0.06 

0.06 

0.06 

0.06 

■m  vvv  rj,--.  vv  iui  --.v-  .mt 


for  the  lift  from  the  mid-section  to  the  tip-section  was 
broadened.  As  tor  the  two  best  models,  the  lift 
distribution  of  model  19  grew  just  very  little  faster  than 
that  of  model  2.  It  also  shows  a  slightly  larger  value  in 
the  tip  region  than  for  the  DLM. 

Again,  the  bumps  in  the  transonic  experimental  pressure 
distributions  were  not  modeled  by  the  numerics.  Also,  the 
discrepancies  between  the  experimental  and  the  numerical 
results  increased  with  increased  Mach  number  and  decreased 
reduced  frequency  due  to  the  linearity  of  the  code,  since 
the  numerical  results  do  not  reflect  shock  waves  possibly 
present  in  the  transonic  cases  M=0. 8  and  0.9. 

The  six  models  were  also  analyzed  in  the  supersonic 
region.  The  lift  coefficient  distributions  are  given  in 
Figures  2b  and  27,  and  in  Tables  2i  through  2b  . 

At  S=1 .1 ,  model  20  shows  the  lowest  lift  from  root  to 
the  80%  span  station  but  shows  a  significant  increase  in 
lift  in  the  tip  region.  Model  2  exhibits  an  almost  linear 
distribution  along  the  span  and  has  less  lift  in  the  tip 
region  than  any  other  model.  Model  19,  which  showed 
reasonably  good  distributions  in  the  subsonic  cases 
experiences  more  lift  than  the  other  endplate  models.  The 
flat  plate  model  shows  less  lift  from  the  root-section  to 
the  aid-section  of  the  wing  but  abruptly  increased  lift  in 
the  tip  region.  The  cruciform  model  and  the  combination 
model  show  high  lift  distributions  over  the  entire  wing. 


TAELE  23 


Unsteady  Litt  Coefficient  Distributions  (Real,  M  =  1.1, 
k=0.11b,  KFM,  6  Models) 


Ming 
Spa  n 

0.0 

(1) 

26.2 

54.1 

75. b 

91.0 

99.0 

Model 

1 

1.27 

1. 52 

1.76 

1.99 

2.19 

2.  38 

Model 

2 

1.26 

1.50 

1.73 

1.92 

2.04 

2.  13 

Model 

4 

1.33 

1. 59 

1.85 

2.10 

2.34 

2.  57 

Model 

14 

1.65 

2.07 

2.74 

3.77 

5.58 

7.73 

Model 

19 

1.  12 

1.57 

1.80 

1.98 

2.09 

2.  1b 

Model 

20 

1.20 

1.45 

1.69 

1.91 

2.12 

2.  33 

TABLE  24 

Unsteady  Litt 

Coefficient  Distributions 
k=0.H6f  KFM,  6  Models) 

(Imag  ,  M=1 .  1 , 

Bing 

0.0  28.2  54.1  75.6 

91.0  99.0 

Span 

<*> 

Model 

1 

0.29 

0.29 

0.29 

0.  29 

0.29 

0.  29 

Model 

2 

0.28 

0.28 

0.28 

0.28 

0.28 

0.  28 

Model 

4 

0.31 

0.31 

0.31 

0.  31 

0.31 

0.31 

Model 

14 

0.55 

0.55 

0.55 

0.  55 

0.55 

0.  55 

Mode  1 

19 

0.22 

0.22 

0.22 

0.  22 

0.22 

0.22 

Mode  1 

20 

0.35 

0.  35 

0.35 

0.  35 

0.35 

0.  35 

Fij'ice  26  : 


'Instead/  Lift  Coefficient  Distributions 
K  =  0.  H6.  KFM ,  6  Models) 


M=1.  1 


TABLE  25 


Unsteady  LiJtt.  Coetticient  Distributions  'Real,  11  =  1.35, 
k=0. 101,  KFM,  6  Models) 


8  ing 
Spa  n 

0.0 

ro 

28.  2 

54.1 

76.6 

91.0 

99.0 

Model 

1 

1.02 

1.19 

1.41 

1.70 

1.94 

2.  19 

Model 

2 

1 . 02 

1.19 

1.40 

1.65 

1.82 

1.98 

Model 

4 

1. 03 

1.21 

1.44 

1.75 

2.04 

2.  32 

Model 

14 

1.04 

1.  1  8 

1.50 

2.  33 

3.76 

5.  47 

Model 

19 

1.02 

1.19 

1.41 

1.67 

1 .86 

2.01 

Model 

20 

1.02 

1.19 

1.40 

1 .67 

1.89 

2.  13 

TABLE  26 

Unsteady  Lift  Coefficient  Distributions  (imag,  M=1.35, 
k=0. 101,  KFM,  6  Models) 


Wing 

Span 

0.0 

:*) 

28.2 

U" 

-P 

• 

-4 

75.6 

91.0 

99.0 

Model 

1 

0.06 

0.06 

0.06 

0.  0b 

0.0b 

0.0b 

Model 

2 

0.06 

0.06 

0.06 

0.06 

0.06 

0.0b 

Model 

4 

0.06 

0.06 

0.0b 

0.  06 

0.06 

0.  06 

Model 

14 

0.07 

0.07 

0.07 

0.  07 

0.0  7 

0.07 

Model 

19 

0.06 

0.06 

0.06 

0.0b 

0.06 

0.06 

Model 

20 

0.06 

0.06 

0.06 

0.  06 

0.06 

0.06 

MODEL  1 
MODEL  2 
MODEL  4 
MODEL  14 
MODEL  19 
MODEL  20 


. . . . I . i 1 1  . . . . . . . 

0.2  0.3  O.li  0.5  0.5  0.7  0.5  0.9  1.0 

y/s 


Unsteady  Lift  Coefficient  Distributions  'M=i.  3o 
K=0.  U  1 ,  KF 'I »  6  Models) 


When  the  flow  became  more  supersonic,  N=1.3b,  the  utt 
distributions  ot  ail  six  models  dropped  by  a  consideraoie 
amount  (20%  at  the  root).  Otherwise,  they  showed  almost  tut 
same  distributions  troa  the  root  section  to  the  mid-section 
and  slightly  increased  lift  in  the  tip  region.  The 
cruciform  still  shows  an  abrupt  increase  starting  at  the 
mid-section.  The  endplate  models,  again,  gave  lower  lift 


distributions  than  the  other  models, 
results  were  available  for  comparison. 


ho  experimental 


The  computational  efforts  for  these  models  by  the  KFM 


are  listed  in  Tables  27  through  32 


As  can  be  seen,  the 


KF3  is  faster  than  the  DLM  and  needs  less  storage.  And  as 
in  the  DLM,  the  endplate  models  take  less  computer  time  and 
storage  than  the  other  models. 


4.4  COHPABISON 

From  sections  4.2  and  4.3,  we  can  deduce  the 
superiority  of  the  cylindrical  slender  body  model  and  the 
endplate  models  over  the  other  models.  Therefore,  in  this 
section,  a  comparison  was  conducted  of  the  endplate  models  2 
and  19,  the  slender  body  model  22,  and  the  flat  plate  model 
1  with  the  results  of  Dusto’s  complete  model  (8)  and  with 
the  experimental  results  of  NLR  (9). 

Here,  the  unsteady  pressure  distributions  in  chordwise 
direction  were  compared  at  M=0.  8  for  sections  2,  b,  and  8. 
These  pressures  are  listed  in  Tables  33  through  J8  and  shown 
in  Figures  28  through  33  . 
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TABLE  27 

Comparison  of  Effort  (M=0.6,  k=0.200,  KFM,  6  Models) 


CPU  Time 
(sec) 

Storage  I/O 
(K-min) 

Disk  I/O 
( EXCP) 

Model 

1 

37.27 

647.5 

1716 

Model 

2 

38.46 

687.1 

1812 

Mode  1 

4 

37.30 

647.9 

1716 

Model 

14 

40.  96 

780.2 

2051 

Mode  1 

19 

37.92 

653.3 

17  15 

Model 

20 

37.60 

650.4 

1715 

Comparison 

of 

TABLE 

Effort  ( M  =  0.  7 

28 

,  k=0. 173 , 

KFM ,  6  Models 

CPU  Time 

Storage  I/O 

Disk  I/O 

(sec) 

(K-ain) 

(EXCP) 

Model 

1 

37.24 

649.7 

1729 

Model 

2 

38.32 

688.5 

1825 

Model 

4 

37.  21 

649.7 

1729 

Model 

14 

40.78 

782.2 

2064 

Model 

19 

37.44 

651.3 

17  28 

Model 

20 

37.  39 

650.8 

1728 

TABLE  31 

Comparison  ot  Effort  (M=1.10,  k=0.116,  KFM,  6  Models) 


CPU  Time 
(sec) 

Storage  I/O 
(K-min) 

Disk  I/O 
( 2XCP) 

Model 

1 

38.5b 

662.3 

1729 

Mode  1 

2 

41.64 

720.2 

1825 

Model 

4 

38.77 

664.3 

1729 

Mode  1 

14 

47.23 

844.5 

2064 

Model 

19 

39.47 

670.8 

1728 

Mode  1 

20 

39.47 

670.9 

17  28 

TABLE  32 

Comparison  of  Effort  (M=1.35,  k=0.l0i,  KFM,  b  Models) 


CPU  Time  Storage  I/O  Disk  I/O 
(sec)  (K-min)  (EXCP) 


Model 

1 

38.14 

655.6 

1717 

Mode  1 

2 

41.16 

711.8 

1812 

Model 

4 

38.21 

656.1 

1717 

Model 

14 

4b. 26 

829.8 

2052 

Model 

19 

38.91 

6  b  2. 4 

1716 

Model 

20 

39.07 

6  6  3.9 

171b 

tfl  vC 
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TABLE  33 

Pressure  Distributions  at  Section  2  (Real,  M  =  0.8,  k=0.1b3) 

for  Various  Models 


TABLE  34 

Pressure  Distributions  at  Section  2  (Iaag,  M=U.3,  !c=D.l53) 

tor  Various  Models 


DLM  1 
DLM  2 
DLM  19 
KFM  1 
KFM  2 
KFM  19 
DUSTO  UP 
DUSTO  LO 
EXP.  UP 
EXP.  LO 


^.DUSTO 


. . I . . . I'" . . . "I . ’"I 

I  0.2  0.3  0.«  0.S  0.0  0.7  0.0  I 


x/c 


Pressure  Distributions  at  Section  2 
*=0.3,  k  =  0.  1  53)  for  Various  'lo-lels 
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TABLE  3b 

Pressure  Distributions  at  Section  5  (Real,  M=0.B,  k.=0.  153) 

tor  Various  Models 


Chord 

Model  1 

Model  2 

Modei19 

Model22 

D  usto 

E  xp 

(*> 

DLM 

KFM 

DLM 

KFM 

DLM 

KFM 

DLM 

«m4 

3 

11.  84 

10. b9 

11.70 

10.43 

11. 68 

10.67 

11.62 

24.32 

0.66 

• 

10 

b.  47 

b.b4 

b.  38 

b.  40 

5.  37 

b.bO 

5.32 

6.04 

4.75 

20 

3.  83 

3.93 

3.7  b 

3.  78 

3.  73 

3.82 

3.67 

3.  54 

3.98 

*  J 

30 

2.88 

3.04 

2.80 

2.90 

2.  80 

2.90 

2.7  3 

2.59 

3.27 

40 

2.  21 

2.37 

2.  14 

2.24 

2.  13 

2.23 

2.06 

2.03 

2.28 

bO 

1.69 

1.82 

1.63 

1.71 

1.  62 

1.68 

1.56 

1.56 

1  .74 

•V 

60 

1. 27 

1.3b 

1.22 

1.25 

1.  22 

1.22 

1.16 

1.14 

1.20 

70 

0.  34 

0.94 

0.89 

0.87 

0.  90 

0.84 

0.8b 

0.82 

0.75 

80 

0.  66 

0  .  b9 

0.62 

0.  54 

0.63 

0.52 

0.59 

0.56 

0.44 

90 

0.41 

0.3  1 

0.38 

0.28 

0.  39 

0.27 

0.36 

0.33 

0.14 

99 

0.22 

0.14 

0.21 

0.01 

0.21 

0.01 

0.20 

0.  1  1 

- 

TABLE  36 

Pressure  Distributions  at  Section  5  (Imag,  M=0.  8,  k=0.i53) 

tor  Various  Models 


Chord 

t*) 

Model  1 
DLM  KFM 

Model  2 
DLM  KFM 

Modell 9 
DLM  KFM 

Model22 

DLM 

Dusto 

Exp 

3 

-1.73 

-2.  69 

-1.65 

-2.  58 

-1.65 

-  1.96 

-1.59 

-6.63 

-0.18 

10 

-0.35 

-1. 03 

-0.31 

-0.  97 

-0.  31 

-0.46 

-0.28 

-1.12 

-0.88 

20 

0.18 

-0.28 

0.  21 

-0.  25 

0.20 

0.22 

0.23 

-0.09 

-0.17 

■* 

30 

0.50 

0.20 

0.51 

0.  21 

0.  5  1 

0.62 

0.52 

0.36 

0.07 

40 

0.70 

0.  53 

0.70 

0.  54 

0.70 

0.88 

0.71 

1.24 

0  .30 

50 

0.82 

0. 7b 

0.82 

0.  74 

0.8  1 

1.01 

0.82 

0.76 

0  .61 

60 

0.86 

0.  85 

0.86 

0.  84 

0.85 

1.04 

0.86 

0.35 

0.b4 

\  • 

70 

0.85 

0.  84 

0.84 

0.  83 

0.  8  4 

0.97 

0.8  4 

0.35 

0.54 

-  *  ] 

80 

0.77 

0.  72 

0.7b 

0.  71 

0.  76 

0.80 

0.75 

0.77 

0.3  7 

90 

0.60 

0.49 

0.59 

0.  49 

0.59 

0.54 

0.58 

0.57 

0.18 

99 

0.38 

0.03 

0.  38 

0.  03 

0.38 

0.03 

0.37 

0.24 

- 

Lti 


TABLE  38 


Pressure  Distributions  at  Section  8  (Imag, 

tor  Various  Models 


n=o.a,  k=o.iS3) 
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In  the  DLM,  since  the  pressures  are  calculated  as  tae 
average  values  of  the  upper  and  lover  surface  pressures,  the 
distribution  ot  pressure  should  be  mbetween  the  two 
experimental  results.  In  the  KF8,  the  pressure 
distributions  are  those  of  the  upper  surface,  so  the  results 
of  the  KFM  are  expected  to  follow  the  shape  of  the  upper 
surface  results  of  the  experimental  data.  As  can  be  seen  Ll 
the  figures,  both  K Fa  and  DLS  show  good  agreement  with  the 
results  of  Dusto,  and  also  show  fairly  good  agreement  with 
experimental  results. 

At  section  2,  the  results  of  the  KF8  show  better 
agreement  than  those  of  the  D Lti  with  the  experimental 
results  except  at  101  -  20%  chord.  At  section  8,  both 
methods  do  not  predict  a  small  pressure  rise  near  the 
trailing  edge.  For  all  cases,  the  real  parts  of  the 
pressure  show  better  comparison  than  the  imaginary  parts. 

As  can  be  seen  in  the  figures,  the  flat  plate  model  and 
the  endplate  models  have  very  similar  pressure  distributions 
except  at  the  tip-section.  In  general,  the  KFfl  displays 
somewhat  more  accurate  distributions  that  the  DLM. 


Chapter  V 

COBCLOSIOS  AMD  POTOBE  BOBK 


5.1  COMCLDSIOH 

Calculations  ot  the  aerodynamic  characteristics  of  an 
F-b  wing  with  and  without  a  wing  tip  missile  were  carried 
out  by  two  computer  codes,  a  doublet  lattice  code  JH7WC)  and 
a  kernel  function  code  (ANKF)  .  The  unsteady  lift 
distributions  in  spanwise  direction  and  the  unsteady 
pressure  distributions  in  chordwise  direction  were  computed 
for  the  wing  with  tip  mounted  launcher  and  missile  without 
fins. 

The  following  conclusions  can  be  drawn  from  the 
performed  analyses: 

1.  In  almost  all  cases,  the  tip  store  can  most 
economically  and  most  accurately  be  represented 
by  an  endplate  or  a  cylindrical  slender  body. 

2.  The  tip  store  increases  the  load  on  the  wing  by 

a  considerable  amount  especially  in  the  tip  region. 

3.  As  the  free  stream  velocity  increases,  the  load  is 
increased  over  the  entire  wing,  but  is  decreased 
when  the  flow  becomes  supersonic. 

4.  The  lift  increases  consistently  with  the  complexity 


of  the  store  model.  Also,  complex  models  require 
more  computer  time  and  storage. 

5.  The  kernel  function  method  is  very  efficient  with 
respect  to  computer  cost.  It  seems  to  be  nearly 
3  times  faster  than  the  doublet  lattice  method 
and  requires  only  about  1/3  of  the  effort  of  the 
doublet  lattice  method  for  all  models.  However, 
the  doublet  lattice  method  seems  to  yield  slightly 
better  results. 

6.  The  endplate  model  2  gives  the  best  results  for 
supersonic  flow,  while  the  endplate  model  19  and 
the  slender  body  model  22  seem  to  be  better  for 
subsonic  flow. 

7.  All  endplate  models  (2,  19,  and  20)  seem  to  be 
better  than  the  present  flat  plate  model  for 
computing  lift  distributions  of  wings  with  tip  stores. 

8.  Finally,  the  two  methods  show  fairly  good  agreement 
with  experimental  results. 

5-2  POTOBB  BOBK 

In  this  paper,  transonic  flow  was  not  treated  as  such. 
It  was  regarded  as  either  a  subsonic  or  supersonic  flow. 
Therefore,  further  research  is  needed  to  account  for 
transonic  flow  effects  such  as  the  shock  wave.  Secondly, 
the  results  from  supersonic  flow  were  not  compared  to 
experiment.  In  the  future,  a  comparison  is  required  to 


ensure  the  validity  of  the  results  for  the  endplate  models. 
Thirdly,  it  was  assumed  that  the  flach  number  distribution  be 
uniform  over  the  entire  wing.  It  is  believed  that  this 
assumption  caused  errors  in  lift  and  pressure  distributions, 
especially  in  the  region  from  M=0.8  to  M=1.1  where  we  can 
expect  some  shock  wave  near  the  trailing  edge.  However, 
this  assumption  was  considered  useful  and  necessary  for  the 
process  of  flutter  computation  when  there  are  no 
experimental  data  for  specific  local  Mach  numbers  over  the 
wing  available. 

Also,  the  differences  in  pressure  distributions  on  the 
upper  and  lower  surfaces  were  neglected.  Since  the  wing  is 
not  symmetric  about  the  camber  line,  the  aerodynamic 
characteristics  of  the  upper  surface  and  lower  surface  of 
the  wing  cannot  be  the  same.  Bore  sophisticated  codes  are 
required  to  account  for  these  differences. 

Finally,  the  aspect  most  important  for  flutter  analysis 
is  the  dynamic  effect  of  the  tip  stores  on  the 
characteristics  of  the  entire  wing.  Thus,  further  research 
is  required  in  the  way  of  performing  flutter  analysis  by  use 
of  large  aeroelastic  codes  such  as  FASTOP  or  SSC/NASTHAN . 
These  analyses  should  compare  the  flutter  results  of  the 
presently  used  flat  plate  models  to  those  obtained  with  the 
aerodynamics  for  the  endplate  and  maybe  slender  body  models 
for  both  KFB  and  DLM.  Additionally,  flight  test  results  or 
wind  tunnel  test  results  should  be  used  for  comparison  and 


evaluation.  If  the  new  nodels  should  prove  advantageous 
similar  analyses  to  the  one  presented  here  could  b 


performed  for  other  tip  stores  and  even  underwing  stores. 

I 
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